function  Direct_Fields( a , b )

GeneRegulation_Nullcline ( a , b );

[u1,v1]=meshgrid(linspace(0,4));
hold on;
streamslice(u1,v1,-u1+a./(1+v1.^b),-v1+a./(1+u1.^b));

end

